Few-Body Systems 0, 1-6 (2008) f^W 

Systems 



© by Springer- Vcrlag 2008 
Printed in Austria 



Bright vortex solitons in Bose Condensates 



Sadhan K. Adhikari* 



Instituto de Fisica Teorica, Universidade Estadual Paulista, 
01405-900 Sao Paulo, Sao Paulo, Brazil 



Abstract. We suggest the possibility of observing and studying bright vortex 
solitons in attractive Bose-Einstein condensates in three dimensions with a ra- 
dial trap. Such systems lie on the verge of critical stability and we discuss the 
conditions of their stability. We study the interaction between two such soli- 
tons. Unlike the text-book solitons in one dimension, the interaction between 
two radially trapped and axially free three-dimensional solitons is inelastic 
in nature and involves exchange of particles and deformation in shape. The 
interaction remains repulsive for all phase 5 between them except for 5 ~ 0. 



1 Introduction 



Solitary waves or solitons are a consequence of nonlinear dynamics. A classic 
text-book example of soliton appears in the following one-dimensional nonlinear 
free Schrodinger equation in dimensionless units 



dt dy 2 



"*737 " 733 " \*M 



V(y,t) = 0. (1) 



The solitons of this equation are localized solution due to the attractive nonlinear 
interaction — \^(y, t)\ with wave function at time t and position y: \P(y,t) = 
\Y2\n\ exp(— ifit)sech(yy/\f2\), with Q the energy pQ. 

Solitons have been noted in optics [2], high-energy physics and water 
waves P|, and more recently in Bose-Einstein condensates (BEC) |U p|. The 
Schrodinger equation with a nonlinear interaction — |<PJ 2 does not sustain a lo- 
calized solitonic solution in three dimensions. However, a radially trapped and 
axially free version of this equation in three dimensions does sustain such a 
bright solitonic solution |Hj which has been observed experimentally Here 
we study the dynamics of these bright solitons. We also suggest that such soli- 
tons can be generated in an axially rotating nonzero angular momentum state, 
which are called bright vortex solitons and can be observed in BEC. 



* E-mail address: adhikari@ift.unesp.br 



2 



Bright vortex solitons in Bose-Einstein Condensates 



A number of bright solitons constituting a soliton train was observed in an 
experiment by Strecker et al. [3], where they turned a repulsive BEC of 7 Li 
atoms attractive by manipulating the background magnetic field near a Feshbach 
resonance [7,. It was found [3] that solitons in such a train usually stay apart. 
Also, often a soliton was found to be missing from a train [3]. There have been 
theoretical attempts [HHH1E3 to simulate essentials of these experiments j3J|3]. 

We use the explicit numerical solution of the axially-symmetric mean-field 
Gross-Pitaevskii (GP) equation ^IJ to study the dynamics of bright solitons in 
a soliton train A j. Attractive BEC's may not form vortices in a thermodynam- 
ically stable state. However, due to the conservation of angular momentum, a 
vortex soliton train could be generated by suddenly changing the inter-atomic 
interaction in an axially-symmetric rotating vortex condensate from repulsive to 
attractive near a Feshbach resonance [T in the same fashion as in the experi- 
ment by Strecker et al. [3J for a non-rotating BEC. Alternatively, a single vortex 
soliton could be prepared and studied in the laboratory by forming a vortex in 
a small repulsive condensate and then making the interaction attractive via a 
Feshbach resonance and subsequently reducing the axial trap slowly. 



2 Mean-field Model and Results 



Mean-field Model: In a quantized vortex state ^2], with each atom having an- 
gular momentum Lh along the axial y axis, the axially-symmetric wave function 
can be written as ^(r, r) = ip(r,y,r) exp(iL6) where 8 is the azimuthal angle 
and r the radial direction. The dynamics of the BEC in an axially-symmetric 
trap can be described by the following GP equation I12j 



d d 



1 d 
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dt dr 2 r dr dy 2 4 ^ ^ ^ ' 
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<p{r,y,t) = 0, 



(2) 



where the length, time, and wave function are expressed in units of \fhj (2muj), 
and [r\J Z 3 /v / 8]~ 1 , respectively. Here radial and axial trap frequencies are 

oj and Xlo, respectively, m is the atomic mass, I = yh/Jmu) is the harmonic 
oscillator length, and n = Na/l the nonlinearity with a the interatomic scat- 
tering length. For solitonic states n is negative. In terms of the one-dimensional 
probability P(y,t) defined by 



P(y,t) = 2tt f dr\tp(r,y,t)\ 2 /r, 
Jo 



(3) 



the normalization of the wave function is given by dyP(y, t) = 1 

We solve the GP equation numerically using split-step time-iteration 
method using the Crank-Nicholson discretization scheme described recently jl.Sj . 
The details of the numerical scheme for this problem can be found in |14j . 

Results: For bright solitons the nonlinearity n is negative. Under the con- 
ditions n < and A = 0, a soliton-type BEC state can be generated only for 
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Figure 1. Three-dimensional wave function \ifi(r, y)/r\ vs. r and y for a single soliton with 
A = 0, and (a) L = 0, n = -0.2, and (b) L = 1, n = -1. 



n greater than a critical value (n C r): n cr < n < 0. For n < n cr , the system 
becomes too attractive and collapses and no stable soliton could be generated. 
The actual value of ncr is a function of the trap parameter A. For the spherically 
symmetric case A = 1, and n C r = —0.575 [TT1 IT2] . 

Numerically solving the GP equation © for A = 0, we find that the critical 
n for collapse of a single soliton is n C r = —0.67 for L = and ncr = —2.10 for 
L = 1. For L = n C r = —0.67 is in close agreement with with ncr = —0.676 
obtained by other workers [HJEil- For L = 1, ncr = —2.10 should be contrasted 
with n C r = —2.20 obrained by Salasnich [T^]. For L = 0, uo = 27rx800 Hz and final 
scattering length — 3ao as in the experiment of Strecker et al. [3], ncr = —0.67 
corresponds to about 6000 7 Li atoms. One can have proportionately about three 
times more atoms in the L = 1 state. 

A L = soliton with n = —0.2 is illustrated in Fig. 1 (a) where we plot the 
three-dimensional wave function \ip(r,y)/r\ vs. r and y. For L = 1 we calculated 
the soliton for n = — 1 and plot \(p(r,y)/r\ in Fig. 1 (b). Because of the radial 
trap the soliton remains confined in the radial direction r, although free to move 
in the axial y direction. The nature of the two wave functions are different. For 
L = 0, the condensate has maximum density for r = 0. For L = 1, because of 
rotation a vortex has been generated along the axial direction corresponding to 
a zero density for r = 0. 

Next we consider the interaction between two solitons with a phase difference 
5 given by the following superposition of two solitons (p at ±yo at time t = 
<f(r,y) = \if(r,y + yo)\+e lS \(p(r,y — yo)\. The time evolution of these two solitons 
is found using the solution of © for different 5. In the present simulation we 
consider two L = 1 equal vortex solitons each of n = —0.4 for A = initially at 
positions yo = ±15 and observe them for an interval of time t = 400. We also 
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Figure 2. The wave function \ip{r,y)/r\ vs. r and y of two L — solitons with n = —0.2 each 
and S = tt/2 at times (a) t = 0, (b) 100, (c) 200 and (d) 300. 



consider the evolution of two L = solitons each of n = —0.2 for yo = ±15. 
The solitons interact by exchanging particles and after an interval of time two 
unequal solitons are generated from two equal solitons. The evolution of the two 
solitons in the L = case for 5 = tt/2 is shown in Fig. 2. Similar evolution for 
the L = 1 case is reported in [Tl] . 

From Fig. 2 we find that at t = the two solitons are equal and symmetri- 
cally located. However, this symmetry is broken for t > 0. The asymmetry and 
separation in the final position of the solitons y\ and 1/2 are best studied via 
I C 1 2/1 1 — 1 2/2 1 ) I and \y\ — 2/2 1 , respectively, at time t = 400 for different phase S 
between the solitons and in Fig. 3 (a) we plot the same for different 5 for the 
L = case. We find that the asymmetry is zero for 5 = tt and and is largest for 
a 5 in between. However the separation increases monotonically as 5 increases 
from to tt. Hence the interaction is repulsive for almost all 5 except for 5 ~ 0. 
Closely associated with the asymmetry and separation in the final position of 
the solitons is the number of exchanged atoms between the two solitons, which 
demonstrates the change in the sizes of the solitons. Actually, the smaller soliton 
travels faster and the larger one travels slower. This results in the asymmetry in 
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Figure 3. (a) Final asymmetry and separation in the position of the two solitons of L = and 
n = —0.2 each at t = 400 vs. phase 8. (b) Ratio Nr/N vs. t for different phase 8 for two L — 
and n — —0.2 solitons initially at ±15. Nr is the number of atoms in the right soliton and N 
the total number of atoms, (c) Same as (a) for two solitons of L = 1 and n = —0.4. (d) Same 
as (b) for two solitons of L = 1 and n = —0.4. 



the final positions. The change in the sizes of the solitons for L = is demon- 
strated in the plot of Nr/N vs. t for different 5 in Fig. 3 (b), where Nr is the 
number of atoms in the right soliton and iV the total number of atoms in the two 
solitons. In Figs. 3 (c) and (d) we plot the same for L = 1. The variation of Nr/N 
is qualitatively similar for L = and 1 in Figs. 3 and also to that found in 
for L = 0. However, there are quantitative differences, specially at large times. 
In the present simulation we find that, for the change 5 — * —5, Nr/N — * Nl/N, 
where Nl = (N — Nr) is the number of atoms in the left soliton. 

If the phase difference 5 between two neighboring solitons is not close to 
zero, they experience overall repulsion and stay apart. However, for 5 close to 
zero they interact attractively and often a soliton could be lost as observed in the 
experiment of Strecker et al. [1]. Throughout this investigation in the interaction 
of two equal solitons we assumed that the nonlinearity \n\ for each is less than 
|ncr|/2, so that a stable solitonic condensate with total \n\ < |ncr| exists when 
the two coalesce. However, if two solitons each with |n| > |ncr|/2 encounter for 
(5 = 0, the system is expected to coalesce, collapse and emit atoms via three- 
body recombination. It is possible that in this case only a smaller single soliton 
survives. This might also explain some missing soliton(s) in experiment. 
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3 Conclusion 

We emphasize the possibility of creating and studying bright vortex solitons 
of an attractive BEC in laboratory under radial trapping. We determine the 
condition of critical stability of bright solitons. Employing a numerical solution 
of the GP equation with axial symmetry, we have performed a realistic mean-field 
study of the interaction among two bright solitons in a train and find the overall 
interaction to be repulsive except for phase 5 between neighbors close to 0. There 
is an inelastic exchange of atoms between two solitons resulting in a change of 
size and shape. Except in the 5 ~ case, the solitons in a train stay apart and 
never cross each other as observed in the experiment by Strecker et al. @]. For 
5 ~ a single soliton can often disappear as a result of the attractive interaction 
among solitons, as observed experimentally by Strecker et al.. The L = 1 vortex 
solitons can accommodate a larger number of atoms and the present study may 
motivate future experiments with them. 
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